***REPLCIATION DO FILE FOR Jessica Fortin, 2013, IPSR "Measuring presidential powers: some pitfalls of agreggate measurement"
***JFR 2014-09-26

*TABLES USING SCand FHT data.dta, FROM SHUGART AND CAREY AND FRYE, HELLMANN and TUCKER
****TABLE 1

***note: one case is different and some figures do not replicate to the second decimal (minimal differences with published paper)
pwcorr paveto plveto decree exintro budget referenda cabinet cabinet2 censure diss total if duplicates==0 & democracy==1, sig

***TABLE 2
factor paveto plveto decree exintro budget referenda cabinet cabinet2 censure diss if  duplicates==0 & democracy==1, mineigen(1)

***diagnostics on these (not in the paper)
scree
***parallel analysis (you have to install package to run paran!!)
paran paveto plveto decree exintro budget referenda cabinet cabinet2 censure diss if duplicates~=1 & democracy==1, iterations(5000)

***TABLE 7
***error in paper, cabinet dismissal and formation are inverted in order, I did not notice it in the proofs, but this does not change the results
factor paveto plveto exintro cabinet cabinet2 censure if duplicates~=1 & democracy==1, mineigen(1)
scree
*** KMO
estat kmo
***Cronbach's Alpha
alpha paveto plveto exintro cabinet cabinet2 censure if duplicates~=1 & democracy==1


***NOW USING SECOND DATA FILE FROM K. ARMINGEON AND CAREJA
****TABLE 3: Patterns of Association between 29 Unweighted Constitutional Powers in 28 Post-communist Countries
***using Tau-b for pairwise, symmetrical associations (ordinal)
ktau ac_1 ac_2 ac_3 ac_4 ac_5 ac_6 ac_7 ac_8 ac_9 ac_10 ac_11 ac_12 ac_13 ac_14 ac_15 ac_16 ac_17 ac_18 ac_19 ac_20 ac_21 ac_22 ac_23 ac_24 ac_25 ac_26 ac_27 ac_28 ac_29 , stats (taub p)

***TABLE 4 Patterns of Association between 29 Weighted Constitutional Powers in 28 Post-communist Countries
***using spearman's rank correlation coefficients
spearman WAC_1 WAC_2 WAC_3 WAC_4 WAC_5 WAC_6 WAC_7 WAC_8 WAC_9 WAC_10 WAC_11 WAC_12 WAC_13 WAC_14 WAC_15 WAC_16 WAC_17 WAC_18 WAC_19 WAC_20 WAC_21 WAC_22 WAC_23 WAC_24 WAC_25 WAC_26 WAC_27 WAC_28 WAC_29 , stats (rho p)

****FIGURE 1: Kernel Density Plot of Unweighted Presidential Power Index
kdensity total, normal


****FIGURE 2: Kernel Density Plot of Weighted Presidential Power Index
kdensity total2, normal


****TABLE 6: Factor Analyses of Unweighted and Weighted Presidential Power Indicators in 28 PostcommunistCountries
paran ac_2 ac_5 ac_6 ac_7 ac_8 ac_9 ac_10 ac_11 ac_13 ac_16 ac_18 ac_19 ac_23 ac_28, iterations(5000)
factor ac_2 ac_5 ac_6 ac_7 ac_8 ac_9 ac_10 ac_11 ac_13 ac_16 ac_18 ac_19 ac_23 ac_28, factors(1)
estat kmo

paran WAC_1 WAC_3 WAC_5 WAC_7 WAC_8 WAC_9 WAC_10 WAC_11 WAC_12 WAC_15 WAC_16 WAC_17 WAC_18 WAC_19 WAC_22 WAC_23 WAC_24 WAC_25 WAC_26 WAC_27 WAC_28 WAC_29 , iterations(5000)
factor WAC_1 WAC_3 WAC_5 WAC_7 WAC_8 WAC_9 WAC_10 WAC_11 WAC_15 WAC_16 WAC_17 WAC_18 WAC_19 WAC_22 WAC_23 WAC_24 WAC_25 WAC_26 WAC_27 WAC_28 WAC_29, factors(2)
rotate, blank(0.4)


***TABLE 8 Factor Loadings of 11 Unweighted Presidential Power Items
***procedure to eliminate via Kaiser-Meyer-Olkin measure of sampling adequacy, repeted iteratively
factor ac_1 ac_2 ac_3 ac_4 ac_5 ac_6 ac_7 ac_8 ac_9 ac_10 ac_11 ac_12 ac_13 ac_14 ac_15 ac_16 ac_17 ac_18 ac_19 ac_20 ac_21 ac_22 ac_23 ac_24 ac_25 ac_26 ac_27 ac_28 ac_29
estat kmo

****.... continue to remove weaker items until final model is selected.
****final model selected
factor ac_5 ac_7 ac_8 ac_9 ac_10 ac_11 ac_16 ac_18 ac_19 ac_23 ac_28, factor(1)
paran ac_5 ac_7 ac_8 ac_9 ac_10 ac_11 ac_16 ac_18 ac_19 ac_23 ac_28, iterations(5000)
estat kmo
alpha ac_5 ac_7 ac_8 ac_9 ac_10 ac_11 ac_16 ac_18 ac_19 ac_23 ac_28

****END, JFR.
